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We discuss the relation between the concentration of the Brown-York mass and the formation of 
trapped surfaces in non-spherical massive systems. In particular, we formulate and prove a precise 
version of the Thorne hoop conjecture in conformally flat three-geometries sliced by equipotential 
foliation leaves. An intriguing relationship between the total rest mass and the Brown-York mass is 
shown. This is a further investigation of the previous work on the Brown-York mass hoop conjecture 
in spherical symmetry. 
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I. INTRODUCTION 

The Brown-York quasi-local mass [l] is one of many 
mass concepts developed in last decades. Its potential 
has been revealed in Q, where the Brown-York mass 
has been used in order to prove the trapped surface 
conjecture [3j in spherically symmetric geometries. The 
trapped surface conjecture states that large mass en¬ 
closed in a small volume has to be trapped, and it con¬ 
stitutes an attempt to concretize a loose idea expressed 
by Thorne in his hoop conjecture Q. 

In this paper we demonstrate that the Brown-York 
mass is useful in proving the hoop conjecture in certain 
classes of non-spherical geometries. They include systems 
having an equipotential surface foliation, that is convex 
in a certain sense. We present sufficient conditions for 
the existence of trapped surfaces. It appears convenient 
to split the consideration into two stages. In the first step 
one deals with 2-surfaces that satisfy an integral condi¬ 
tion - that they are averaged trapped surfaces [s]. In the 
second step one finds additional conditions that ensure 
the pointwise trapping of averaged trapped surfaces. 

The paper is organized as follows. Section [TT] gives 
a concise historical account on the hoop conjecture. 
Section mu contains the description of the formalism 
and needed definitions. We propose various sufficient 
and necessary conditions for the existence of averaged 
trapped surfaces in general settings. This is done in terms 
of the reference geometry. Section IIVI defines an equipo¬ 
tential foliation, assuming conformally flat geometry. We 
provide therein necessary and sufficient conditions for av¬ 
eraged trapped surfaces in terms of the original physical 
geometry. As an aside, but important result, we show 
that the Brown-York mass is not larger than the total 
rest mass. The question whether an averaged trapped 
surface is indeed pointwise trapped is examined in Sec- 
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tion El There is an extra mass term required to balance 
the non-sphericality, that can be written down explicitly. 
The last Section summarizes obtained results. 


II. HOOP CONJECTURE 

There is a folk belief in general relativity that if mat¬ 
ter is sufficiently concentrated into a finite volume, the 
gravitational system ultimately has to collapse to a black 
hole. Thorne proposed a hoop conjecture (HC) [1] which 
states: 

Horizons form when and only when a mass 
M gets compacted into a region whose cir¬ 
cumference in EVERY direction satisfies C < 

At: Ad. 

His conjecture deals with global event horizons and the 
‘circumference’ and the mass are deliberately left unspec¬ 
ified. Notice, however, that in the Schwarzschild space- 
time we have the equality: C = An Ad. In this case M is 
the asymptotic mass [6]. Seifert Q formulated the more 
concrete trapped surface conjecture (TSC), according to 
which massive singularities have to be surrounded by a 
closed trapped 2-surface. This is an easier concept, be¬ 
cause trapped surfaces are local in time. Proving the HC 
would require the study of the full history of a space- 
time, while in order to prove the TSC one needs only to 
consider a single Cauchy slice. 

There have been many attempts to prove the HC/TSC. 
In the early period the concentration of matter was as¬ 
sumed in spherically symmetric spacetimes Early 

results have been reviewed in [lfj. Recently Khuri (l7l | 
applied in this context the generalized Jang equation 
[l8| |. Schoen and Yau dealt with nonsymmetric space- 
times 0- Their sufficient condition for the formation 
of trapped surfaces required a special spacelike foliation 
of a spacetime with large extrinsic curvature. Within a 
single Cauchy slice, assuming the matter density to be 
large on a ’large region’, trapped surfaces have to form. 
That was a consequence of the blow-up analysis of the 
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Jan g eq uation [20] for an asymptotically flat initial data 
set |2l|. In later studies the TSC has been proved in 
special classes of systems, with matter 0, El, IH or in 
vacuum [U,[24j|. Some of the recent development has been 
reviewed in [251 ]. 

In spherical symmetric systems one can measure their 
’size’ by the circumference. It is reasonable to take 
C = 2nR where R is the Schwarzschild or the areal ra¬ 
dius of the surface in question, i.e. R = y 7 Area/47r. 
Then one can prove a precise statement of the hoop con¬ 
jecture using the Brown-York mass as the mass measure 
[3]. The theorem says that if C < 27nriBY) then the sur¬ 
face is trapped. There exists a discrepancy between the 
47 t in Thorne’s HC and the coefficient 2-7 t in Q- This 
can be traced back to the fact that at the horizon of the 
Schwarzschild spacetime the Brown-York mass is equal 
to R, the area radius of the horizon. That is utby = 2 M, 
where M is the asymptotic mass. 


III. AVERAGED TRAPPED SURFACES 

Let (fl 3 , g, K) be a subset of a Cauchy slice for the Ein¬ 
stein field equations. Here g is the 3-metric of a Cauchy 
hypersurface and K stands for its extrinsic curvature. We 
assume that H is time-symmetric, i.e. it lies in a totally 
geodesic Cauchy hypersurface, K = 0. 

In this paper, we concentrate on the case of being 
a compact domain with boundary and since it is time 
symmetric, we use the Brown-York mass as our measure 
of the mass within Cl. 

Assume further that the boundary E = dCl is a topo¬ 
logical 2-sphere. There exists a unit normal n (directed 
outward) to E; its divergence V^n 1 is equal to the mean 
curvature k. Here V, denotes the covariant derivative 
with respect to the 3-metric g. The sign of k has an 
important physical meaning. Take a bundle of outgoing 
null rays, normal to E. If k > 0 along E, then the bundle 
is divergent; while if k < 0, then the null rays must con¬ 
verge. If k < 0 everywhere along E, then the two-surface 
is said to be trapped. Trapped surfaces do not exist in 
the Euclidean geometry and their presence is associated 
with strongly curved geometries. 

On the other hand, the first derivative of the area of 
E with respect to the uniform normal deformation gives 
the total mean curvature, 

H( E) = J fcdE. (1) 

The concept of a trapped surface is purely local, but it 
appears useful to deal with surfaces that are trapped in 
the average: 

Definition A surface E is called an averaged trapped sur¬ 
face (ATS) if ff(E) is negative. 

Assume further that E has positive Gauss curvature 
and thus can be isometrically embedded into the Eu¬ 
clidean space R 3 , i.e. * : E i(E) C R 3 . This isometric 


embedding is called the Weyl embedding and it is unique 
up to a rigid motion in R 3 [261 ]. 

Then the Brown-York mass [l] is defined as 

™by(£,s) = ^ J (fc 0 - fc)dE (2) 

where k is the mean curvature of E with respect to the 
physical metric g and k 0 is that of *(E) with respect to the 
Euclidean metric. Note that ko is completely determined 
by the intrinsic 2-metric on the surface E but does not 
depend on the extrinsic geometry how E bends in Cl. 

The above definition implies in a straightforward way 
the important proposition. 

Proposition 1 The surface S is an ATS if and only if 

Wby(E ,g) > J fc 0 dE. (3) 

It is interesting that here the integral and the mean cur¬ 
vature k 0 in Proposition 1 are in the Euclidean space. 
One can employ well known geometric estimates and re¬ 
express the proposition in a number of ways. This is 
done in the remainder of this section. The total mean 
curvature f s fc 0 dE represents an ’averaged size’ of a solid 
convex body in the Euclidean space [27|. Suppose that 
a compact oriented convex surface E lies in M 3 . Let xq 
be a fixed point enclosed by E. The Minkowski integral 
formula [281, Lemma 6.2.9, Page 136] gives 

/>E = J K(x) < n(x),X(x) — Xq >r 3 dE. (4) 

Here K(x) is the Gauss curvature and X(x) is the posi¬ 
tion vector of E in R 3 , n(x) is the unit normal at X(x) 
and < •, • > denotes the Euclidean inner product. 

Recall that E is a topological sphere. By the Gauss- 
Bonnet theorem f y K(x)d£ = 2wx(S 2 ) = 4n, it 
gives an upper bound of the right hand side of 0, 
47rsup s6S \X(x) -x 0 \. 

If we measure the ’size’ of a surface by looking at the 
position vector of its image when embedded isometrically 
into R 3 , then we have 

Theorem 1 (Sufficient Condition for an ATS.) If 

m BY (E, 5 ) > sup \X(x) - x 0 |, (5) 

xGS 

then E is an ATS. 

Another upper bound of the total mean curvature is 
given by the Blaschke cap body inequality, cf. Page 387 
in [29]. Let V be a compact convex body in R 3 . Then 


Area(dU) > ^ 3Vol(V) J yd(<9U) 

where ko is the mean curvature of the boundary dV. This 
leads to the following theorem. 
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Theorem 2 (Sufficient Condition for an ATS.) If 


m B Y(£,s) > 


1 (Area(E)) 2 

47t 3Vol(flo) 1 


( 6 ) 


then E is an ATS. Here flo is the body in R 3 enclosed by 
the image ofT, via the (unique) Weyl embedding. 

At this stage we have a hybrid picture. The Brown- 
York mass lives in a physical space while the upper 
bounds are given in the reference space. Again things 
are easy in the spherically symmetric case 0, when 
ko = 2 /R, where R is the areal radius. In general Rie- 
mannian geometries, life becomes harder. It is difficult 
to define a workable concept of a ’circumference ’ OS, ( 33 . 
Fortunately, one finds a quantitative link between the 
total mean curvature in the reference space k odE and 
the original physical geometric data in a class of folia¬ 
tions of conformally flat 3-manifolds. The details will be 
discussed in the next section. 

To provide necessary conditions for an ATS, we need 
the lower bound estimates of the total mean curvature 
f s fcodE. There are two candidates both of which are in 
terms of intrinsic 2-geometry of the surface. One is given 
by the classical geometric inequality [29] and the other 
one is give n by the Birkhoff invariant of the intrinsic 2- 
metric [32l |. 


Theorem 3 (Necessary Condition for an ATS.) Assume 
that E is an ATS, then 


(1) m BY (E(7) 

(2) m B y(E,s) > -4/3 = (8) 

Here (3 is the Birkhoff invariant of the surface E. 


IV. CONFORMALLY FLAT GEOMETRIES 

Herein we shall investigate the following concrete class 
of three-spaces. Assume that 

(1) g is conformally flat, g ab = f 4 g ab where g ab is the 
standard Euclidean metric. 

( 2 ) There is an equipotential foliation on ft, 

9 = f A {a)[g a adcr 2 + gijdx'dx^} ( i,j = 2 , 3 ) ( 9 ) 

where a > 0 and a foliates the level surfaces of / which 
are assumed to be convex, and x 2 and x 3 are quasi-angle 
variables. 

( 3 ) E = {a = a 0 }. 

Thus, n = (h CT , 0,0) and h a = y/g a <r- The conformal 
factor / satisfies the elliptic equation A/ = —2npf 5 , 
which is the Hamiltonian constraint for momentarily 


static initial data of the Einstein equations. The energy 
density p is nonnegative due to energy conditions. 
Remark: In order to detect whether a surface E is 
trapped or not, one only needs the geometry in a neigh¬ 
borhood of E in 12, or within E. The dominant view 
nowadays is that trapped surfaces are of physical inter¬ 
est because of their roles in the proof of scenarios of the 
cosmic censorship. That demands that H constitutes a 
domain of an asymptotically flat Cauchy slice and hence 
one should assume that the equipotential surface folia¬ 
tion, that covers SI, is extendible onto the entire slice 
with the asymptotic condition /(oo) = 1. That in turn 
implies that /|s > 1 by the maximum principle. 

We emphasize that E refers to the a = oo surface with 
induced metric f 4 (ao)(gijdx l dx 3 ). Denote by E the a- 
constant surface with induced metric gijdx z dx J . Let k 
be the mean curvature of E with respect to the physi¬ 
cal metric g and let k be the mean curvature of E with 
respect to the Euclidean metric g acr da 2 + gijdx'dxf. 

We isometrically embed E into the reference space 
f 4 (ao)[g aa dcr 2 +gijdx l dx 3 ] which is also Euclidean. Then 
it gives a relation between ko and k , i.e. ko = k/ / 2 (oo) 
and that of the induced area forms is dE = / 4 (oo)dE. 

Below we shall write down some criteria for ATS’s, 
obtained in Section uni in terms of the geometry of the 
physical space (S 71, g). 

(i) Proposition [Tj states that E must be an ATS if 
m B y(E, g) > f s &odE; but f E fc 0 dE = / 2 (oo) f E k dE rep¬ 
resents an ’averaged areal size’ Rav of a body enclosed 
by the (convex) 2-surface E. The sufficiency condition 
states simply m B y(E,g) > Ra v (E). 

(ii) In the same way one may also rewrite Theo¬ 
rem □ as: If m BY (E ,g) > R sup ( E), where R sup ( E) := 
/ 2 (cto) sup^-gf; |A'(x) — ieo|, then E is an ATS. Here X 
is the position vector of the surface E in the Euclidean 
space g aa do 2 + g^jdx’dx 3 . 

As a consequence of the uniqueness of the Weyl em¬ 
bedding, the lower bounds given in Theorems [5] and [3] are 
completely determined by the intrinsic 2-geometry on the 
surface. They are the same no matter calculated either 
in the physical space (f2, g) or in the reference Euclidean 
space. In particular, 

(iii) If we define the areal radius of E as Rs := 
y / Area(E)/47r, then the first condition in Theorem [3] be¬ 
comes m B y(E,g) > R$. 

(iv) In the second condition of Theorem [3] the Birkhoff 
invariant /? is the minimum len gth of a closed string 
being slipped over the 2-surface [33(. One defines the 
Birkhoff radius Rb := /3/27r and then the condition be¬ 
comes rriBY(E,g) > Rg. 

We have introduced the above four size measures for E. 
Notice that in spherical geometries and for a round sphere 
E centered at the symmetry center, all these measures 
coincide, Ra v = Rsup = Rs = Rb. 

The physical and ’embedded’ mean curvatures along E 






4 


are related as 


k\t(a- 


412 =-Pi!? + 7 < I0 > 


There is a simple calculation that allows us to give an 
upper bound onto the Brown-York mass triBY(kff) by 
the total rest mass M(E) = f Q pdvol g within E. Indeed, 


mBY(£, 9 ) = T^_i_ /Vo)d £ 

-r,IjpW) + pk)’ TV ' lf)fMdE 

= [ f( a o)pf 5 f~ 6 dvol g = [ p^^-dvolg 
Jn Jn J 

( 11 ) 

where A/ = —2irpf 5 . If we assume the dominant (or 
weak) energy condition, i.e. p > 0, then / is a superhar¬ 
monic function (with respect to the Euclidean metric) 
and the maximum principle yields that for any x € 17, 
f(x) > /|s = /(oo). Therefore, 

Theorem 4 One has meY(S, g) < M(E) where M(E) = 
In pdvolg. 

Remark: The total rest mass has been employed in 

[3, 1, H3, HH in the derivation of sufficient conditions for 
ATS’s and further TS’s under certain additional condi¬ 
tions. The conditions therein are of the form M( E) > 
D( E) where D(E) is a certain ’size measure’ coming 
from an upper bound of the geometric size of the do¬ 
main enclosed by E. As a corollary of Theorem [H if 
rriBY(£!,<?) > -D(E), then E must be an ATS or TS. One 
is expecting to find a refined size measure D'(Yf (which 
is smaller than D(E)) for the Brown-York mass. We shall 
do it in the next section. 

If E is a marginally trapped massive shell where the 
derivative of / has a discontinuity, then Eq. ED , as a 
consequence of integration by parts, is no longer valid. 
But in spherical symmetry, one shows that the total rest 
mass equals twice of the asymptotic mass [35!, Eq.(1.14)]. 
This value also agrees with the Brown-York mass. 


V. FROM AVERAGED TRAPPED SURFACE 
TO TRAPPED SURFACE 

In spherically symmetric geometries, if we take a spher¬ 
ical two-surface E centered at the symmetry center, then 
its mean curvature becomes a constant. That means that 
it is trapped if and only if it is an ATS. This is the situ¬ 
ation considered in |2|. In non-spherical geometries, the 
existence of an ATS is not sufficient to make use of the 
Penrose singularity theorem [36| to conclude that a black 
hole has to develop. In this section, we formulate cer¬ 
tain additional conditions which guarantee that an ’el¬ 
lipsoidal’ ATS is indeed trapped. We would ask how 


much Brown-York mass compacted into the system can 
produce a pointwise TS. The reasoning is analogous to 
that used in jToj, [H, [34]. 

Suppose that E is not a TS, there must be at least 
a point on which the mean curvature k is nonnegative. 
Then the maximal value of n a k must be nonnegative and 
hence 


J n a (n a k) max dT, > 0. (12) 

There is no summation for a here. Instead, n a denotes 
the particular er-component of the unit normal in the 
, equipotential foliation ©• Then one must have 

E J n a [(n a k) max -n <7 k]dT, + -^ J fc 0 dE > m B Y(£,s)- 

(13) 

Equivalently, we have 

Proposition 2 If 

m B Y(£,sO > J n a [(n a k) rnax -n C rk]d'E + -^-- J k 0 dE, 

^ E (14) 

then E must be a pointwise TS. 

Now we apply Eq.(12) in [23]: 


d a ( / fcdE) = 2 / Kh a dE := 8irC(a). (15) 
Js it 


Integrating from 0 to ctq, we have 


1 


1 


^U k ^ = - 


8tt Jt(a=a 0 ) Pi 0 o) 


/V o)dE 


/ 2 M 

8tt 


/E(<r=0) 


kE-t-/Vo) / C(s)ds. 

(16) 


Note that the {a = 0} ’surface’ is the set of points for 
which the conformal factor / achieves its maximal value. 

One may find an upper bound of the first term, Eq.(18) 
in (37|, 


/ 2 M 

87 r 


JfcdE< 


/ 2 (0-0)7 


/S(<7=0) 


supZ(5(0)). (17) 


Here sup /(;S'(0)) is the largest flat radius of the disk on 
which the conformal factor / achieves its maximal value. 
Finally, we have arrived at 

Theorem 5 (Sufficient condition for a pointwise TS.) If 
m by(£,sO > fff) j na i( n ° k ) max Tier k] dE 


/ 2 (ct o)tt 


supZ(S(0)) 


( 18 ) 


r<* 0 

+ / 2 (ct 0 ) / C(s)ds, 
then E is a pointwise TS. 
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The physical significance of this theorem is as follows. 

(1) The line integral term / 2 (cr 0 ) f^° C(s)ds represents 
an appropriate ’size’ of the hoop for mass concentration. 

(2) The surface integral term / s ^[(nu^moi — 
rio-fc]dE/87r reflects the ’boundary effect’ when the sur¬ 
face is not spherical. 

(3) The radius term f 2 (ao)n sup Z(5(0))/4 is influenced 
by the behavior of the conformal factor / within the en¬ 
tire foliation and thus can be interpreted as the ’global 
deviation’ the system from being spherical. 

The sum 

^ J n°[(n a k)max-n a k\d'E+ f sup 1(5(0)) (19) 

is the energy required to balance the non-sphericality 
when producing a trapped surface. By maximum prin¬ 
ciple, for 0 < a < <7o> /(c) > /(<t o). Then the 
radius term f 2 (ao)ir sup 2(5(0))/4 and the hoop term 
/ 2 (ct 0 ) /q CT ° C(s)ds are both less than the terms in the 
sufficient condition for TS in terms of M(E) [23], 
/ 2 (0)7rsupi(S(0))/4 := 7rrad(0)/4 and f°° f 2 (s)C(s)ds, 
respectively. However, the ’boundary effect’ energy com¬ 
pensating terms J- S n a [(n (7 k) max — n a k]dY^/8'n are the 
same. We obtain a sufficient condition for TS’s, employ¬ 
ing the Brown-York mass, that is finer than that implied 
by Theorem @] (cf. Remark beneath Theorem []]) ■ 


In spherical symmetry, both of the two terms in Eq. 
flUJ) vanish, and C(s) = 1 and hence / 2 (er 0 ) J^° C(s)ds 
equals the areal radius of the surface. The inequality (fTSl) 
is sharp and it reduces to the result in 

VI. CONCLUSIONS 

We have shown, using conformally flat geometries and 
a suitable foliation, that the Brown-York mass is bounded 
from above by the total rest mass. The more non- 
spherical is a surface, the more Brown-York energy must 
be compacted within to make it trapped. We employ a 
number of geometric inequalities in Euclidean space, that 
yield several necessary and sufficient conditions for ATS’s 
and pointwise trapped surfaces. These results hold true 
for a large class of non-spherical geometries whose met¬ 
rics are conformal (with convex layer surfaces) to the flat 
metric, and for adapted foliations. 


Acknowledgments 

N. Xie is partially supported by the National Sci¬ 
ence Foundation of China (grants 11171328, 11121101, 
11421061) and CSC Program. 


[1] J.D. Brown and J.W. York, Phys. Rev. D 47, 1407 
(1993). 

[2] N. O Murchadha, R.-S. Tung, N. Xie, and E. Malec, 
Phys. Rev. Lett. 104, 041101 (2010). 

[3] H.J. Seifert, Gen. Relativ. Grav. 10, 1065 (1979). 

[4] K.S. Thorne in Magic without Magic ed. .J. Klauder (Free¬ 
man, San Francisco), 231-258 (1972). 

[5] J. Hartle and D. Wilkins, Phys. Rev. Lett. 31, 60 (1973). 

[6] R. Arnowitt, S. Deser, and C.W. Misner, ‘The dynamics 
of General Relativity’ in Gravitation: an introduction to 
current research , Ed. L. Witten, (Wiley, NY, 1962); Gen. 
Relativ. Grav. 40, 1997 (2008). 

[7] P. Bizon, E. Malec, and N. 0 Murchadha, Phys. Rev. 
Lett. 61, 1147 (1988). 

[8] P. Bizon, E. Malec, and N. O Murchadha, Class. Quan¬ 
tum Grav. 6, 961 (1989). 

[9] D. Eardley, J. Math. Phys. 36, 3004 (1995). 

[10] J. Guven and N. O Murchadha, Phys. Rev. D 56, 7658 
(1997). 

[11] J. Guven and N. O Murchadha, Phys. Rev. D 56, 7666 
(1997). 

[12] E. Malec, Phys. Rev. D 49, 6475 (1994). 

[13] E. Malec and O Murchadha, Class. Quantum Grav. 21, 
5777 (2004). 

[14] T. Zannias, Phys. Rev. D 45, 2998 (1992). 

[15] T. Zannias, Phys. Rev. D 47, 1448 (1993). 

[16] E. Malec, Acta Phys. Pol. B 22, 829 (1991). 

[17] M. Khuri, Phys. Rev. D 80, 124025 (2009). 

[18] H. Bray and M. Khuri, Asian J. Math. 15, 557 (2011). 

[19] R. Schoen and S.-T. Yau, Commun. Math. Phys. 90, 575 


(1983). 

[20] P.S. Jang, J. Math. Phys. 19, 1152 (1978). 

[21] R. Schoen, and S.-T. Yau, Commun. Math. Phys. 79, 
231 (1981). 

[22] P. Koc and E. Malec, Acta Phys. Pol. B 23, 123 (1992). 

[23] E. Malec, Phys. Rev. Lett. 67, 949 (1991). 

[24] R. Beig and N. O Murchadha, Phys. Rev. Lett. 66, 2421 
(1991). 

[25] L. Szabados, Living Rev. Relativity 12, 4 (2009). 

[26] L. Nirenberg, Comm. Pure Appl. Math. 6, 337 (1953). 

[27] G. Polya and G. Szego, Isoperimetric inequalities in 
mathematical physics. Annals of Mathematics Studies 27, 
(Princeton University Press, Princeton, 1951). 

[28] W. Klingenberg, A course in differential geometry, 
Translated from the German by David Hoffman. Grad¬ 
uate Texts in Mathematics 51, (Springer-Verlag, New 
York-Heidelberg, 1978). 

[29] R. Schneider, Convex Bodies: The Brunn-Minkowski 
Theory, (Cambridge University Press, Cambridge, 2013). 

[30] E. Flanagan, Phys. Rev. D 44, 2409 (1991). 

[31] J.M.M. Senovilla, Europhys. Lett. 81, 20004 (2008). 

[32] J.C.A. Paiva,Bull. Belg. Math. Soc. Simon Stevin 4, 373 
(1997). 

[33] G.D. Birkhoff, Trans. Amer. Math. Soc. 18, 199 (1917). 

[34] P. Bizon and E. Malec, Phys. Rev. D 40, 2559 (1989). 

[35] P. Bizon, E. Malec, and N. O Murchadha, Class. Quan¬ 
tum Grav. 7, 1953 (1990). 

[36] R. Penrose, Phys. Rev. Lett. 14, 57 (1965). 

[37] E. Malec, Acta Phys. Pol. B 22, 347 (1991). 




